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Different from the two-dimensional (2D) topological insulator, the 2D topological crystalline insulator (TCI) 
phase disappears when the mirror symmetry is broken, e.g., upon placing on a substrate. Here, based on a new 
family of 2D TCIs — SnTe and PbTe monolayers — we theoretically predict the realization of the quantum 
anomalous Hall effect with Chern number C — 2 even when the mirror symmetry is broken. Remarkably, 
we also demonstrate that the considered materials retain their large-gap topological properties in quantum well 
structures obtained by sandwiching the monolayers between NaCl layers. Our results demonstrate that the TCIs 
can serve as a seed for observing robust topologically non-trivial phases. 

PACS numbers: 73.22.-f, 73.21.Fg, 73.43.-f 


The topological crystalline insulator (TCI) is a recently dis¬ 
covered new class of quantum states with insulating bulk en¬ 
ergy gap and gapless edge or surface states Q. Gapless edge 
or surface states of TCIs arise from the crystalline symmetry 
in contrast to that of the Z 2 topological insulator (TI), where 
the reason for protection is time-reversal symmetry (TRS) [Q]- 
4]. The first proposal of a material class with TCIs was 
purely theoretically 0 and consequently TCIs were experi¬ 
mentally observed in three-dimensional (3D) IV-VI semicon¬ 
ductors SnTe 0, Pbi-^Sn^Se 0, and Pb^Si^Te CUD. 
Later on, several transition metal oxides imiHl, Bi 2 Te 3 CD, 
and heavy-fermion compounds lfl2l[l3l have been found to be 
3D TCIs. 

Rather recently, the TCI state has been theoretically pre¬ 
dicted to occur in two dimensions (2D) in multilayer struc¬ 
tures of SnTe fPfl . PbSe rm. and graphene multilayers GE 
While the experimental progress in 3D TCIs was rapid m- 
Eiin), the 2D TCI phase disappears when placed on a sub¬ 
strate for further experimental investigation or device appli¬ 
cation due to the breaking of the mirror symmetry. Without 
the substrate, the freestanding film is usually hard to grow, 
which is also the main reason why 2D TIs were experimen¬ 
tally mainly established in HgTe/CdTe and InAs/GaSb quan¬ 
tum wells up to now mum. Therefore, it is highly desirable 
to search for new materials and novel formation mechanisms 
which maintain the mirror symmetry of 2D TCIs. 

TCIs bear great potential for the investigation of exotic phe¬ 
nomena Ii l 7 , 2014221 . such as the large Chern number quantum 
anomalous Hall (QAH) effect 121,22|. The QAH state, which 
is characterized by a quantized charge conductivity without an 
external magnetic field, has topologically protected dissipa¬ 
tionless chiral edge states and therefore is a good starting point 
to realise ultra-low-power consumption electronics |231. The 
TRS breaking in TIs provides a promising platform to investi¬ 
gate this striking topological phenomenon (241125) . However, 
the Chern number of the Tl-based QAH state is usually lim¬ 
ited to C = 1 1 24-261 or C = 2 (271 128) . In magnetically 
doped thin-film TCI SnTe, the QAH phase was predicted with 
the Chern number reaching as much as four El 122) . The 
QAH phase with large Chern number leads to multiple dis¬ 


sipationless edge channels that provide better ways to opti¬ 
mise electrical transport properties and could greatly reduce 
the contact resistance for circuit interconnects 12111281 291. 

Here, based on density functional theory, we predict that 
SnTe and PbTe (001) monolayers are 2D TCIs. We expose our 
TCIs to an exchange field, which mimics the effect of mag¬ 
netic adatoms or a substrate: While an in-plane exchange field 
destroys the TCI phase, it survives for an out of plane magne¬ 
tization. Based on the considered materials, we demonstrate 
the feasibility of realizing a TCI-originated QAH phase with 
high Chern number, which is very robust to mirror symme¬ 
try breaking and/or magnetization fluctuations. Remarkably, a 
range of 2D TCIs can be realised in NaCl 2 Ar(Sn/Pb)Tei quan¬ 
tum wells with a single layer of SnTe or PbTe embedded in 
N layers of widely used substrate NaCl i3Ql on both sides. 
These systems exhibit large topologically non-trivial energy 
gaps with a magnitude tunable by the NaCl thickness. Our 
findings pave the way to utilization of TCIs as seed materials 
for arriving at robust topologically non-trivial phases. 

The first-principles calculations have been performed using 
the full-potential linearized augmented-plane-wave method as 
implemented in the FLEUR code ED- The self-consistent 
calculations with spin-orbit coupling (SOC) were carried out 
with a cutoff parameter, fc max , of 3.8 bohr -1 . The experi¬ 
mental lattice parameters were used for pristine Sn/PbTe sys¬ 
tems, while those for quantum well structures were obtained 
using the Vienna ab-initio simulation package (VASP) (32). 
The generalized gradient approximation in parametrization of 
Perdew, Burke and Ernzerhof (GGA-PBE) was used for the 
exchange correlation potential (33). The maximally local¬ 
ized Wannier functions (MLWFs) were constructed using the 
wannier90 code Hi- 

Bulk SnTe and PbTe share the same face-centered cubic 
NaCl-type structure with an inverted band structure at the L 
point, which results in the realisation of a 3D TCI state in 
bulk SnTe and a 2D TCI state in SnTe multilayers (T4). On 
the other hand, PbTe both in bulk and in thin (001) films is in 
the normal insulator (NI) state |2 # [14) . Here, we focus on the 
(001) oriented monolayers of SnTe and PbTe. In Fig.[lja) the 
top view of such 2D monolayer is shown, with Sn (Pb) and Te 
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FIG. 1. (color online) (a) Top view of the crystal structure of (001) 
Sn/PbTe monolayers with the unit cell indicated by dashed lines, (b) 
Schematic view of the NaCl 2 jv(Sn/Pb)Tei quantum well structure 
for N = 2. (c) Corresponding 2D and projected ID Brillouin zones 
with marked high symmetry points, (d) Berry curvature distribution 
associated with d =i mirror eigenstates of the occupied bands along 
the T — X — T path. 


atoms forming two square sublattices positioned in the mirror 
plane z = 0. Accordingly, all Bloch states in the system can 
be labeled with the eigenvalues d =i of the reflection operator 
with respect to this symmetry plane. 

To get a preliminary insight into the topological proper¬ 
ties of the systems, we present in Fig. [2] the orbitally-resolved 
band structures of SnTe and PbTe monolayers with and with¬ 
out SOC. In the absence of SOC for SnTe, energy bands with 
the Sn -s and T e-p y orbital character (positive parity with the 
inversion center at the Sn atom) overlap around the X-point 
with the Sn -p z states (negative parity). For PbTe at the Ap¬ 
point without SOC a direct band gap appears with the valence 
band maximum (VBM) and the conduction band minimum 
(CBM) dominated by the Pb -s and T e-p y orbitals (positive 
parity), and the Pb -p z orbital (negative parity), respectively. 
Turning on SOC leads to insulating character in both systems 
(calculated band gaps are 0.05 eV for SnTe and 0.09 eV for 
PbTe), and to the band inversion in PbTe, so that with SOC 
the band structure is inverted in both systems. 

Owing to an even number of X points in the Brillouin 
zone (see Fig.[ljc)) neither the SnTe monolayer nor the PbTe 
monolayer is a 2D TI. However, taking the mirror symmetry 
into account, similarly to bulk SnTe 0, for both SnTe and 
PbTe monolayers band inversion results in the realisation of 
a 2D TCI state. To show this explicitly, we calculate the so- 
called mirror Chern number um = {ji+i ~ n -i)/ 2, where n+i 
and n_i are Chern numbers of all occupied bands with op¬ 
posite mirror eigenvalues +i and — i, respectively f36l . The 
Chern number of a given subset of states is given by C m 
X f BZ Q(k) d 1 2 k , where D(k) is the Berry curvature 0711381 : 


fi ( k ) = 2Im 

u<Ef m^n 


(V^nk \v x | V^mk) (^mk \Vy |^nk) 
(^rak ^nk)^ 


,(D 


with m,nas band indices which run over considered subset of 



FIG. 2. (color online) Orbitally-resolved band structures for (001) 
monolayers of SnTe (a, b) and PbTe (c, d) without (a, c) and with (b, 
d) SOC, weighted with the contribution of s, p y , and p z states. Pari¬ 
ties of the conduction band minimum and valence band maximum at 
the X point are labeled by “+” and The Fermi level is indicated 
with a dashed line. 


states, i/j n k and e n k are corresponding Bloch states and their 
eigenenergies, and v x / y as the cartesian components of the 
velocity operator. In Fig. |TJd), we plot the distribution of the 
Berry curvature of all occupied bands with mirror eigenvalue 
=b i. The main contribution to the Berry curvature comes from 
the region around X , with its values having opposite sign for 
opposite eigenvalues. The Chern number for each polariza¬ 
tion is, respectively, n+i = — 2 and ri-i = 2, yielding the 
total Chern number of all occupied states C = 0 and the mir¬ 
ror Chern number um = —2. The calculated mirror Chern 
number um = — 2 proves the TCI nature of (001) oriented 
SnTe and PbTe monolayers. [] 

One of the prominent features of a TCI is the existence of 
gapless edge states that appear as a result of the crystal mir¬ 
ror symmetry only and do not require TRS mm. To illus¬ 
trate this, using description in terms of MLWFs of 2D mono- 
layers, we construct the tight-binding Hamiltonian of finite 
80-atom wide ribbons oriented along the (llO)-axis. In addi¬ 
tion, by computing the matrix elements of the Pauli matrices 
a a (a = x , y , z) in the basis of MLWFs, we consider the ef¬ 
fect of an exchange field applied to our ribbon, which breaks 
the TRS and generally has out of plane (cr z • B± term in the 
Hamiltonian) and in-plane ( a x ■ B y term) components. By 
introducing such an exchange field we aim at mimicking the 
effect of interaction with a magnetic environment, e.g. mag¬ 
netic adatoms or magnetic substrate l39l . 

The spin-resolved band dispersion of states localized at the 


1 The fact that in Ref. I'f4]| only films thicker than four layers were found to 
be TCIs might be traced back to the fact that the film bandstructures were 
modeled from bulk DFT results in that work. Here, DFT calculations were 

performed explicitly for the films. 
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FIG. 3. (color online) Evolution of the edge states of a (110)-oriented 
ribbon of PbTe monolayer in an applied exchange field, (a) Upper- 
edge without exchange field (B = 0.00 eV, 2D TCI), (b) Upper-edge 
with exchange field within the mirror plane (Bn — 0.08 eV, NI), and 
(c) upper-edge with exchange field perpendicular to the mirror plane 
(B± =0.20 eV, QAHI). Edge states are colored with the expectation 
value of cr z , as shown in (a), (d) Phase diagram of 2D PbTe with 
respect to B\\ and B±. The color scale represents the magnitude of 
the energy gap at the X point (in eV). The dashed line marks the 
phase boundary and the dots correspond to situations shown in (a)- 
(c). The inset in (d) shows the chiral edge states at the upper-edge of 
the ribbon at a point in the phase diagram marked with a dot next to 
it. 


“upper” (“lower”) edge and the projected bulk band struc¬ 
ture without any exchange field is presented in (the inset of) 
Fig-Ha). Clearly, according to the mirror Chern number anal¬ 
ysis, two pairs of edge states cross slightly away from the 
time-reversal invariant X-point, and the spin polarization of 
the crossing states is of the same sign, which is in contrast to 
the case of 2D TIs. For an in-plane exchange field the mir¬ 
ror symmetry is broken and the edge states become gapped 
(Fig.[3](b)). If the exchange field is, on the other hand, out of 
plane, the mirror symmetry survives and the system continues 
to exhibit gapless edge states, although the TRS is broken. 

The TRS breaking causes an exchange splitting between 
conduction and valence bands of opposite spin, and thus, with 
increasing strength of the exchange field they approach each 
other. As shown in Fig. [3jd), with increasing the magnitude 
of B± the band gap decreases with the system remaining a 
TCI, and it closes at around B± = 0.16 eV, reopening again 
with further increasing B±. According to our calculations, 
this signals a phase transition from a 2D TCI to a QAH phase, 
of which the latter is characterized by a non-zero Chern num¬ 
ber of C = +2. The nonzero Chern number is further verified 
by the presence of two gapless edge states on each side of the 
ribbon within the nontrivial band gap (see Fig.[3jc)). 

Our most interesting finding is that the QAH state in our 
TCI system can be reached irrespective of whether the mirror 
symmetry is broken by B\\ or not (see computed phase dia¬ 


FIG. 4. (color online) Orbitally-resolved relativistic band structures 
of (a) NaCUivSnTei and (b) NaCUjvPbTei. Color from red to blue 
represents the increasing contributions of Sn/Pb-s and Tc-p y orbitals, 
which reveals the band inversion clearly, (c) Calculated energy gap 
of NaCl 2 iv(Sn/Pb)Tei superlattice at the X point with and with¬ 
out SOC as a function of NaCl thickness. Negative and positive 
values correspond to the normal and inverted energy gaps, respec¬ 
tively. Inset shows the global energy gap as a function of NaCl thick¬ 
ness. (d) Localization-resolved band structure of a ID nanoribbon of 
NaCUPhTei. Inset is the zoom-in at the X point. 


gram in Fig.^d) and the inset of it). This means that TCIs can 
be used to arrive at novel Chern insulators despite their appar¬ 
ent fragility to the mirror symmetry breaking, e.g., by a sub¬ 
strate. In fact, it is only the magnitude of the exchange field 
and not its direction, which the evolution of the topological 
band gap is sensitive to, except for the situation with B± = 0, 
when the system always resides in the NI state. Such a pecu¬ 
liar behavior also means that, when the QAH state is achieved 
in the system, it is not only robust with respect to the crystal 
symmetry, but also to the direction magnetization of adatoms 
or the substrate, which is subject to thermal fluctuations. 

Finally, we choose the widely used material NaCl lf30l 
to propose a practical way of stabilizing the 2D TCIs while 
keeping the mirror symmetry in NaCl 2 Ar(Sn/Pb)Tei quantum 
wells. One of our main findings here is that the magnitude 
of the TCI band gap can be engineered by the thickness of 
NaCl spacer, and it can reach as much as 0.13 eV and 0.47 
eV in NaCl 2 ArSnTei and NaCl 2 ArPbTei wells, respectively. 
The sketch of the crystal structure of a NaCl 2 Ar(Sn/Pb)Tei 
quantum well is shown in Fig. [T] (b). A SnTe (PbTe) mono- 
layer is placed in between TV-layers thick NaCl films, so that 
the mirror symmetry survives. When SnTe (PbTe) monolayer 
is sandwiched between NaCl slabs, the interaction between 
SnTe (PbTe) and NaCl pushes the Sn/Pb -p z orbital up in en¬ 
ergy and Sn/Pb-s with Te-p y orbitals down in energy at the 
X point. For the case of N = 1 the strongest interaction is 
observed, resulting in the disappearance of the band inversion 
and the TCI phase. With increasing thickness of NaCl layers, 
the strength of the local potential decreases. Consequently, 
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the band inversion and the TCI phase return at a critical NaCl 
thickness. 

To show the thickness-dependent band inversion, we plot 
the orbitally-resolved band structures in Fig. |4|a,b). We ob¬ 
serve that for small N the Sn/Pb-s and T e-p y orbitals con¬ 
tribute to the VBM with positive parity, while, at larger NaCl 
thickness, they contribute to the CBM keeping positive parity. 
Therefore, the band inversion occurs at N = 2 and N = 3 for 
the NaCl 2 ArSnTei and NaCl 2 ArPbTei, respectively. Fig.|4jc) 
shows the magnitude of the energy gaps at the X point without 
and with SOC versus the NaCl thickness. For NaCl 2 ArSnTei, 
the existence of the normal energy gap in both cases demon¬ 
strates a normal insulator phase for TV = 1. For N = 2, 
a SOC-induced band inversion appears, indicating the realisa¬ 
tion of the 2D TCI phase. As the NaCl thickness increases fur¬ 
ther, the energy gap at the X point increases and the inverted 
energy gap is obtained even without SOC. Similar behaviour 
is seen for NaCl 2 ArPbTei with varying N, but the influence 
of SOC is more pronounced in this case. The 2D TCI state 
is further explicitly confirmed by the topological analysis and 
the emergence of the gapless edge states in ID nanoribbons 
(see Fig. |4jd)). The global energy gap is one of the defin¬ 
ing properties of TCI. As shown in the inset of Fig. |4|c), the 
global gap survives for a large range of NaCl thickness, and 
its magnitude can be tuned by a proper choice of N, reach¬ 
ing as much as 0.13 eV and 0.47 eV for NaCl 6 SnTei and 
NaCl 22 PbTei, respectively. We note that the latter value of 
0.47 eV is the largest for either 2D or 3D TCIs reported so far. 

In summary, we have identified by first-principles calcula¬ 
tions that the (001) oriented SnTe and PbTe monolayers are 
topological crystalline insulators. We explore the phase dia¬ 
gram of these materials with respect to an applied exchange 
field due to magnetic adatoms/substrate. We show that as the 
mirror symmetry remains for the out of plane magnetization, 
the gapless edge states survive, but the TCI state is destroyed 
for an in-plane exchange field. For finite B ±, upon reaching a 
certain strength of the exchange field we show the emergence 
of a robust Chem insulator phase which survives irrespective 
of the mirror symmetry breaking. We further show that the 
large-gap TCI Pb/SnTe based family can be successfully sta¬ 
bilised in NaCl 2 Ar(Sn/Pb)Tei quantum wells. Especially for 
NaCl 22 PbTei, an energy gap of 0.47 eV appears, which is by 
far larger than that in known 2D and 3D TCIs. 
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